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Monte Carlo computer simulations are used to study transient cavities and the solvation of hard- 
spheroid solutes in dipolar hard sphere solvents. The probability distribution of spheroidal cavities in 
the solvent is shown to be well described by a Gaussian function, and the variations of fit parameters 
with cavity elongation and solvent properties are analyzed. The excess chemical potentials of hard- 
spheroid solutes with aspect ratios x in the range 1/5 < x < 5, and with volumes between one 
and twenty times that of a solvent molecule, are presented. It is shown that for a given molecular 
volume and solvent dipole moment (or temperature) a spherical solute has the lowest excess chemical 
potential and hence the highest solubility, while a prolate solute with aspect ratio x should be 
more soluble than an oblate solute with aspect ratio 1/x. For a given solute molecule, the excess 
chemical potential increases with increasing temperature; this same trend is observed in the case 
of hydrophobic solvation. To help interpret the simulation results, comparison is made with a 
scaled-particle theory that requires prior knowledge of a solute-solvent interfacial tension and the 
pure-solvent equation of state, which parameters are obtained from simulation results for spherical 
solutes. The theory shows excellent agreement with simulation results over the whole range of solute 
elongations considered. 



I. INTRODUCTION 

The solvation of solutes in simple polar solvents is of 
importance in almost all areas of chemistry and biochem- 
istry. The way in which solvent molecules are ordered 
(or not) around solute molecules can confer all manner 
of interesting and useful phenomena upon the solution, 
including catalytic activity, preferential solvation in so- 
lute mixtures, microphase formation, and self-assembly. 
Some of the most important - and complex - solvation 
phenomena occur in aqueous solutions of hydrophobic so- 
lutes. Hydrophobic solvation is characterized by negative 
energy and entropy changes accompanying the insertion 
of hydrophobic solutes in to aqueous phases. One con- 
tribution to the energy change is from the ever-present 
attractive dispersion interactions that operate between 
all species, but the dominant contribution is from the 
specific solvent-solvent interactions facilitated by a re- 
organization (ordering) of the solvent molecules around 
each solute molecule; this latter effect also gives rise to 
the negative entropy change. As a result of these changes 
in energy and entropy, the solubility of a hydrophobic 
solute in aqueous solution decreases with increasing tem- 
perature over quite a wide temperature range; in the case 
of methane in water the solubility reaches a minimum 
at T ~ 350 K. The hydrophobic attraction between hy- 
drophobic solutes in aqueous solution is an effective inter- 
action arising from the fact that favorable solvent-solvent 
interactions are maximized when the solute molecules are 
in close proximity with one another, as opposed to being 
completely surrounded by solvent. For a recent discus- 
sion of these concepts, and references to the enormous 
literature in this area, see 

Considerable effort has been directed toward under- 
standing the nature of solvation in aqueous solutions, es- 
pecially as compared to that in simple non-aqueous sol- 



vents. For example, inert gases are less soluble in water 
than they are in non-aqueous molecular Iiquidsi2*2i4 So- 
phisticated molecular theories have been proposed specif- 
ically for hydrophobic solvation and attractions^ It was 
pointed out that the distribution of transient cavities 
in the pure solvent could be used to characterize such 
effects as hydrophobic solvationS^ In the early 1990's, 
computer simulations were used to characterize the tran- 
sient cavities in a range of molecular liquids^S*ii This 
was achieved by calculating the probability of finding a 
spherical solute-sized cavity in the solvent. It was found 
that although water has a larger free volume - or al- 
ternatively, a smaller effective packing fraction - than 
do non-aqueous solvents, that free volume is distributed 
throughout smaller cavities. Moreover, non-aqueous sol- 
vents exhibit a greater ability to redistribute free volume 
in order to accommodate large solute species. Such con- 
siderations are of fundamental importance in understand- 
ing the hydrophobic interaction and all of its manifesta- 
tions in biochemistry, nanoscale systems, and materials 
chemistry. 

This very brief discussion highlights the value of un- 
derstanding the molecular-scale cavity structure in the 
pure solvents before one proceeds to consider specific 
solvent-solute systems. In Ref. it was suggested that 
such investigations could be extended to non-spherical 
cavities and solutes. Simulation studies of this type 
have been carried out for various hard-core solute/solvent 
systemsiSiiiiiiiSiiLii but we are not aware of any sys- 
tematic comparison between rod-like and disk-like so- 
lutes. The current work is therefore concerned with 
transient non-spherical cavities and the solvation of non- 
spherical solutes in model polar solvents. For simplicity, 
we consider hard-spheroid solutes (or uniaxial hard el- 
lipsoids) in a solvent made up of dipolar hard spheres 
(DHSs). The DHS system is the simplest model of a po- 
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lar liquid, and is of considerable intrinsic interest owing to 
its complex phase behaviori 18 i 19 i 20 i 21 i 22 i 23 i 24 i 25 i 26 In par- 
ticular, the low-temperature properties of the DHS fluid 
are dominated by the association of particles in 'nose- 
to-tail' conformations giving rise to chains and rings at 
low densities, and extended networks at intermediate 
densitiesiiSiSL 2 ^ Therefore, the DHS system is an inter- 
esting example of a fluid that should exhibit large tran- 
sient cavities as compared to non-polar hard-sphere fluids 
at the same density. We use Monte Carlo (MC) computer 
simulations to study the size distributions of prolate (rod- 
like) and oblate (disk-like) spheroidal cavities in DHS flu- 
ids, and relate these distributions to the excess chemical 
potentials of hard-spheroid solutes. This allows an ex- 
amination of the role of solute shape, as well as solute 
volume, on the solubilities of non-polar solutes in polar 
fluids. Following earlier simulation studies of this type, 
we compare our calculations with a scaled-particle theory 
(SPT) which describes the reversible work for insertion 
and subsequent growth of a single solute molecule in to 
a solvent. 

This article is organized as follows. In Section [H] we 
describe the solute and solvent models, the SPT as ap- 
plied to the model system, and the computer simulation 
methods employed in this work. Results are presented in 
Section ITTT1 and Section llVl concludes the paper. 



II. MODEL AND METHODS 
A. Model 

The DHS solvent consists of hard spheres with diam- 
eter a, each carrying a central dipole moment \x. For 
two spheres with separation vector r, the pair potential 
energy is 
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where r = \r\ and eo is the vacuum dielectric permitiv- 
ity Thermodynamic parameters of the DHS system are 
expressed in dimensionless form as follows: the reduced 
number density p* — Na 3 /V, where N is the number 
of particles in a volume V; the reduced dipole moment 
pf = \J p 2 /AireoksTa 3 , where ks is Boltzmann's con- 
stant, and T is the temperature; the reduced tempera- 
ture T* = {l/p*f. Throughout the rest of the paper, 
the DHSs will be referred to as component '1'. 

We will compute the probability of finding a spheroidal 
cavity in the DHS fluid centered at a randomly selected 
point in the system, and hence obtain the excess chemical 
potentials of hard-spheroid solutes. The solute will be re- 
ferred to as component '2' of the resulting solution. With 
the spheroid radial semi-axes denoted by a, and the polar 
semi-axis by c, the elongation of the spheroid is given by 
x = c/a. The fundamental measures of the spheroid will 
be required in Section [II CI and are therefore recorded in 
Table H] for reference. 



B. Thermodynamics 

Of central interest in this paper is the excess chemi- 
cal potential, A/i 2 , of hard-spheroid solutes in DHS sol- 
vents at low concentrations where solute-solute interac- 
tions can be ignored. Ap 2 can be related to the solu- 
bility of the solute in a variety of situations, but to pro- 
vide a concrete example consider the transfer of a solute 
molecule from an ideal gas to an initially pure solvent 
(or an extremely dilute solution) such that the volumes 
of the gas and the solvent separately remain constant. 
Thermal equilibrium between the gas and the solution 
is assumed. The chemical potential of the solute in the 
ideal gas is equal to 



M2 aP = k 



B T\np v 2 ap V 



(2) 



where p^ p is the number density of solutes in the gas 
phase, and V is the de Broglie thermal volume of the 
solute. The chemical potential of solute molecules in 
solution (in the absence of solute-solute interactions) is 
written 



/i2 = k B T In p 2 V + Ap 2 



(3) 



where p 2 is the number density of solutes in the solution. 
When the solute is at chemical equilibrium the chemical 
potentials of the solute in both phases are equal, and from 
the mass-action law the associated partition coefficient 
can be defined as 
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where f3 — l/kgT. The temperature dependence of the 
partition coefficient yields information on the energetic 
and entropic contributions to (3Ap 2 . For transfer under 
conditions of constant volume the relevant expressions 
are 



Au 2 = k B T 
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for the energy, and 
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(6) 



for the entropy, where all differentiations are carried out 
with gas volume, solvent volume, and number of solvent 
molecules held constant. Widom et al. have provided full 
details of analogous expressions for transfer under condi- 
tions of constant pressure, and the physical situations to 
which Eqs. JSJ and © are applicable iSi 



C. Scaled particle theory 

We note that there are several theories of hard so- 
lute/solvent systems that rely on representing non- 
spherical components by spherical particles with effective 
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TABLE I: Mean radius of curvature, surface area, and volume for prolate and oblate spheroids with radial semi-axis a, polar 
semi-axis c, and elongation x = c/a. 



Prolate (x > 1) 



Eccentricity e 
Mean radius of curvature R2 



a/T~— x~" 



Oblate (a; < 1) 




Surface area 


S 2 


2na 2 


arcsin e 


2na 2 


l+i#ln(i±f)' 


Volume 


v 2 






4 2 
^7ra c 



radii^iSiiSii^iiSiiSiii the resulting effective hard-sphere 
mixture is then described by accurate semi-empirical for- 
mulae, such as the Boublik-Monsoori-Carnahan-Starling- 
Leland equation of statei^SS In this work, we will show 
that the variations of A[i 2 with solute size and shape can 
be described most intuitively within a simple version of 
SPT. 30 The excess chemical potential of a single hard 
spheroid with radial semi-axis a and polar semi-axis c 
is equal to the reversible work required to insert a point 
particle in to the solvent, and then to grow the point par- 
ticle to full size. Consider, then, a scaled spheroid with 
radial and polar semi-axes equal to Xa and Ac, respec- 
tively, with A = corresponding to the point particle, 
and A = 1 to the full-sized particle. For hard-particle 
solutes and solvents, the Widom formula relates the ex- 
cess chemical potential for the solute at reciprocal tem- 
perature [3 = 1/ksT to the probability, P, of inserting 
a solute at random in to the solution without incurring 
any overlaps >2i 
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In the spirit of SPT, we will interpolate between results 
that are correct in the limits A = and A = 00 to obtain a 
result for A = 1. When A is small, the chemical potential 
of the scaled spheroid is related to the free volume in the 
solvent, i.e. 
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where V cx (\) is the excluded volume of the scaled 
spheroid and a single solvent molecule. Obviously, this is 
only correct when the solute is sufficiently small that it 
can only be in the vicinity of one solvent molecule at any 
given time. The sphere-scaled spheroid excluded volume 
is given by the following expression due to Kiharai 3 ^ 



V ex (A) = Vi + XS 1 R 2 + X 2 R 1 S 2 + X 3 V 2 



(9) 



The fundamental measures for spheroids are given in Ta- 
ble [U for DHSs the measures are R\ = a/2, Si = na 2 , 
V± = 7r<7 3 /6. In the other extreme, A — > 00, the reversible 
work for inserting the solute particle is dominated by the 
work done against the pressure, p, of the solvent, and the 
(normally positive) contribution to the free energy aris- 
ing from the interface between the solute molecule and 



the solvent. In this case the excess chemical potential is 

X 2 1 S 2 X 3 P V 2 
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where 7 is the solute-solvent interfacial tension. We now 
approximate the excess chemical potential of a full-sized 
solute molecule with a cubic function of A: the terms 
up to order A are obtained by Taylor expansion of the 
combination of (jHJ and JHJ; and the terms in A 2 and A 3 
are taken from l|10(l . Finally, we set A = 1 (corresponding 
to the full-sized solute molecule) which yields 
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where r\ = pV\ is the packing fraction of the pure sol- 
vent. The first term in Ijlll) arises from the fact that 
the probability of inserting a point particle in to the sol- 
vent without overlap is equal to (1 — tj). (The surface 
term proportional to can also contain an additional 
curvature contribution 3 ^ but it will be explained below 
that this is unnecessary in the current application.) The 
pressure and solute-solvent interfacial tension will be de- 
termined from computer simulations which we describe 
in the next section. 



D. Computer simulations 

Canonical (NVT) MC simulations of N = 500 DHSs 
were carried out in a cubic simulation cell with peri- 
odic boundary conditions applied. 34 The long-range dipo- 
lar interaction was handled using the Ewald summation 
with conducting boundary conditions. One MC cycle 
consisted of one trial translation and one trial rotation 
per molecule, on average. Run lengths after equilibration 
consisted of O(10 5 ) MC cycles. The respective maximum 
displacements for translational and rotational moves were 
adjusted to give acceptance rates of 20%. 

To characterize the cavity structure in the DHS fluid, 
a test particle was inserted in to the simulation cell with 
randomly selected position and orientation. For a given 
elongation, the maximum possible spheroid radial semi- 
axis, Qjfn, that would not result in any particle overlaps 
was determined. This procedure was repeated 500 times 
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every 20 MC cycles, and a histogram of a m was accu- 
mulated to yield a cavity function Q(a m ), the integral of 
which is 



Q(a m )da 

r> 



1 



(12) 



In addition, the excess chemical potentials /3A//2 = 
— In P(a) were determined using J7J) and the fact that 
for a given test-particle insertion, all spheroids with the 
same elongation and a < a m would be accommodated in 
the solvent without overlap. P(a) and Q(a m ) are related 
because the probability of successfully inserting a test 
particle with radial semi-axis a is equal to the likelihood 
of there being a cavity that can accommodate a particle 
with radial semi-axis of at least that sizejiSi i.e. 

/•CfO />a 

P(a)= / Q(a m )da m = 1 - r\ - I Q(a m )da m . (13) 

J a JO 

As will be explained in Section fill HI MC simulation re- 
sults were used to parametrize the SPT expression (|11H . 
The pressure in the DHS fluid was determined using the 
virial equation^ 
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where g(cr) is the radial distribution function, g(r), at 
contact, and U is the dipolar configurational energy. The 
second term in (|14f> represents the hard-core contribution 
to the equation of state; g{a) was estimated by extrapo- 
lating g(r) to contact using a [2,2] Pade approximant. 



III. RESULTS 

The DHS fluid was studied with reduced dipolc mo- 
ments p* = (hard spheres), //* = 1, and p* = 2, cor- 
responding to reduced temperatures of T* = oo, T* = 1, 
and T* = 0.25, respectively. The number densities con- 
sidered were p* — 0.2, p* — 0.5, and p* = 0.8. (For water 
under ambient conditions, p* ~ 2 and p* ~ 0.9; rough es- 
timates of the parameters for liquid ammonia are p* ~ 1 
and p* ~ 0.8.) Hard-spheroid solute particles were con- 
sidered with elongations in the range 1/5 < x < 5, and 
with a range of molecular volumes up to 20 Vi, depend- 
ing on the density; at higher densities the probability of 
inserting a large solute particle in to the solvent without 
overlap becomes too small to enable an accurate evalua- 
tion of J7J). 



A. Cavity distributions and molecular structure 

In Fig. ^ we show the cavity distributions, Q{a m ), for 
spheroidal cavities in a DHS solvent with p* = 0.2 and 
p* = 2. Results are presented for spheroids with elonga- 
tions of x — 1/5 (oblate), x = 1 (spherical), and x = 5 




FIG. 1: Linear-log plot of the cavity distribution functions, 
Q(a m ), for spheroidal cavities of radial semi-axis a and elon- 
gation i in a DHS solvent with p* = 0.2 and p* — 2: x = 1/5 
(circles); x = 1 (squares); x = 5 (diamonds). The solid lines 
are fits to the Gaussian distribution. 



(prolate). The distributions are well described by the 
Gaussian function 



Q{a m ) oc exp 



ao 



(15) 



where ao is the most probable radial semi-axis, and s a is 
the corresponding width parameter. Distributions of the 
maximum polar semi-axis are denned analogously. The 
Gaussian function produced excellent fits to all of the 
DHS state-points and cavity elongations studied; results 
for p* — 0.2 and p* — 2 are included in Fig.^ In Figs.EJ 
13 and0|we show the fitting parameters for DHS solvents 
at densities of p* = 0.2, p* = 0.5, and p* = 0.8, respec- 
tively, as functions of the spheroid elongation. These 
figures also include the corresponding parameters for the 
polar semi-axis, c. 

The results for p* = 0.2 (Fig. |2Jl show that with all 
dipole moments, the most probable radial (polar) semi- 
axis decreases (increases) with increasing elongation. In 
addition, the most probable polar semi-axis for prolate 
cavities with elongation x > 1 is greater than the most 
probable radial semi-axis for oblate cavities with elon- 
gation 1/x < 1; this is just due to fact that the polar 
semi-axis of a long, thin spheroid is longer than the radial 
semi-axis of a short, fat spheroid with the same volume. 
For a given cavity elongation, ao and Co are almost in- 
dependent of the solvent polarity (except at the highest 
elongations considered where small deviations are appar- 
ent). The width parameters in the Gaussian functions 
vary in a similar way to the corresponding most prob- 
able values, although they appear to be more sensitive 
to the solvent dipole moment. Results for p,* = and 
p* = 1 are almost identical, but those for p* = 2 are sig- 
nificantly larger for both semi-axes and all elongations. 
At low dipole moments (or high temperatures) the distri- 
bution of cavity sizes is narrower since there is only one 
characteristic lengthscale in the solvent. At high dipole 
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FIG. 2: Statistics for the radial semi-axis (open symbols) 
and polar semi-axis (filled symbols) of spheroidal cavities in 
DHS solvents at density p* = 0.2 and dipole moments n* — 
(circles), fi* = 1 (squares), and p* — 2 (diamonds): (left) 
most probable values, ao and Co; (right) width parameters, s a 
and s c . Error bars are smaller than the symbols. 



FIG. 3: Statistics for the radial semi-axis (open symbols) 
and polar semi-axis (filled symbols) of spheroidal cavities in 
DHS solvents at density p* = 0.5 and dipole moments /i* = 
(circles), y,* = 1 (squares), and fj,* — 2 (diamonds): (left) 
most probable values, ao and Co; (right) width parameters, s a 
and s c . Error bars are smaller than the symbols. 



moments there is strong association of the solvent parti- 
cles which gives rise to two structural lengthscales: the 
average distance between nearest-neighbor solvent par- 
ticles, and the characteristic lengthscale of the cluster 
network^ This gives rise to a heterogeneous cavity struc- 
ture in which small cavities are located in the vicinity of 
clustered solvent particles, and larger voids are formed in 
the free space between clusters. Hence, a broader distri- 
bution of cavity sizes is found at higher dipole moments. 
The width parameter is clearly a sensitive probe of the 
solvent structure. 

Results for p* = 0.5 (Fig-EJ and p* = 0.8 (Fig.© show 
similar trends, except that the most probable dimensions 
show local maxima in the range of elongations considered. 
For a given solvent dipole moment, these maximal values 
of ao and cq occur for less anisotropic shapes as the den- 
sity is increased, which suggests that there are fewer long 
'holes' in the fluid structure at high densities. The most 
probable cavity dimensions and associated width param- 
eters decrease with increasing density. This is because 
as the density is increased, the free volume in the vicin- 
ity of the clustered solvent molecules represents a greater 
proportion of the total free volume; since small cavities 
are located in the region of the solvent molecules, and 
large cavities are accommodated in the voids of the clus- 
tered network structure, the average cavity dimensions 
and associated 'standard deviations' will be reduced, ao 
and Co decrease with increasing solvent dipole moment, 
whereas s a and s c increase, once again due to the pro- 
motion of solvent-particle association. The development 
of a heterogeneous fluid structure gives rise to a greater 
spread of cavity sizes, whereas the shifts in ao and Cq 
must be due to the decreasing nearest-neighbor separa- 
tion between solvent particles. 

The influences of solvent dipole moment on the fit pa- 
rameters are seen to decrease in magnitude with increas- 




FIG. 4: Statistics for the radial semi-axis (open symbols) 
and polar semi-axis (filled symbols) of spheroidal cavities in 
DHS solvents at density p* = 0.8 and dipole moments /i* = 
(circles), fi* = 1 (squares), and p* — 2 (diamonds): (left) 
most probable values, ao and Co; (right) width parameters, s a 
and s c . Error bars are smaller than the symbols. 



ing density. This is due to the fact that at low densi- 
ties strongly interacting DHSs associate to form networks 
that possess heterogeneous cavity structures, whereas at 
high densities the fluid structure is not qualitatively dif- 
ferent from that of pure hard spheres. This statement 
is backed up by Fig. which shows the static structure 
factor, S(q), in DHS solvents at each density and dipole 
moment considered in this work. There are many alter- 
native choices of structural probe, but 5(g) does have 
the merit of observing certain characteristic scaling laws 
in the presence of strong particle association^ 5(g) was 
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FIG. 5: Log- log plots of the static structure factors, S(q), for 
DHS solvents with dipole moments p* = (circles), /i* = 1 
(squares), and /x* = 2 (diamonds): (a) p* = 0.2; (b) p* = 0.5; 
(c) p* = 0.8. The dashed line in (a) follows the power law 
S(q) OC-T 1 . 



calculated in the simulations directly using the relation 
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(16) 

where q is a reciprocal lattice vector of the simulation 
cell, and contributions with equal q = \q\ were averaged. 
At a density of p* = 0.2 the structure of the system 
with p* = 2 is clearly different from the structures of 
those systems with pi* = I and p* = 0. In particular, the 
strongly polar system exhibits power-law behavior at low 
q; as discussed fully in Ref.|22, chain-like correlations give 
rise to the scaling S(q) ~ q^ 1 . At the higher densities, 
the structural features of the fluid are far less sensitive 
to the dipole moment. 



B. Parametrizing the SPT 

To test the SPT expression in Ijllfl we need to deter- 
mine the values of the solute-solvent interfacial tension 
and the solvent equation of state. The latter was ob- 
tained as described in Section III Dl and the results are 
included in Table |nj The interfacial tension was ob- 
tained by fitting to simulation results for spheri- 




FIG. 6: Excess chemical potentials of hard-sphere solutes 
with radius a in DHS solvents with p* = 2 and at various 
densities: p* = 0.2 (circles); p* = 0.5 (squares); p* = 0.8 
(diamonds). The symbols are from MC simulations and the 
curves are best fits using 111 IB 



cal solutes only, with the intention of then applying the 
parametrized equation to non-spherical solutes. As ex- 
amples of the procedure, in Fig. we show the excess 
chemical potentials of hard-sphere solutes as functions of 
the radius, a, in DHS systems with pf = 2. It is em- 
phasized that the only fitting parameter in (|ll|l is the 
reduced interfacial tension jcr 2 /ksT^ The fits included 
in Fig. are seen to be quite good, and the resulting fit 
parameters arc reported in Table ITT1 In earlier works^S*^ 
it has sometimes proven necessary to modify the interfa- 
cial term with a curvature correction, i.e. the third term 
on the right-hand side of (|11|) is replaced by a contribu- 
tion likeS 



7^2 / AS 
k B T \ R 2 



(17) 



where S is also to be parametrized against simulation re- 
sults. This modification was fully tested against simula- 
tion results, but in some cases it led to negative interfacial 
tensions, and in every case gave rise to significant uncer- 
tainties in the fit parameters with no visible improvement 
over Ullf). 



C. Excess chemical potentials of spheroid solutes 

With i|ll|) parametrized using simulation results for 
spherical solutes, we now turn to an examination of the 
excess chemical potentials for non-spherical solutes. To 
emphasize the effects of molecular shape we compare re- 
sults for hard-spheroid solutes with the same volume, V2, 
but with different elongations (and hence different R2 
and S2). In Fig. 0we present simulation results for DHS 
solvents at p* = 0.2, and solutes with molecular volumes 
in the range V% < V2 < 20 V% . For a given molecular vol- 
ume the spherical solute has the lowest value of f3Ap2, 
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TABLE II: Thermodynamic properties of the DHS fluid: the reduced density p"; the reduced dipole moment p*\ the reduced 
temperature T*; the radial distribution function at contact g(a); the dipolar energy, U, per molecule in units of fcsT"; the 
reduced pressure pa^/ksT; the reduced solute-solvent interfacial tension 7a 2 /fcgT. 



p* 




T* 




U/Nk B T 


pa li /k B T 




0.2 





oo 


1.3220 





0.3107 


0.13103(2) 


0.2 


1 


1 


1.6492 


-0.2847 


0.2812 


0.13085(2) 


0.2 


2 


0.25 


9.6583 


-4.6588 


0.0774 


0.09579(7) 


0.5 





oo 


2.1558 





1.6288 


0.5128(1) 


0.5 


1 


1 


2.4148 


-0.6625 


1.4332 


0.5089(1) 


0.5 


2 


0.25 


5.6678 


-5.7950 


0.5702 


0.4516(1) 


0.8 





oo 


4.0094 





6.1742 


1.3653(3) 


0.8 


1 


1 


4.2377 


-1.0056 


5.6758 


1.3623(6) 


0.8 


2 


0.25 


6.1733 


-6.8986 


3.5558 


1.3153(4) 



and hence from 10} the highest solubility. This reflects 
the well-known fact that the ratio of surface area to vol- 
ume is smallest for spherical molecules, and hence the 
unfavorable interfacial tension term in (|ll|l is minimized 
with this geometry. (3Ap2 increases sharply as the solute 
geometry deviates from the sphere, and there is an ap- 
proximate correspondence between the excess chemical 
potentials of oblate solutes with elongation 1/x < 1 and 
prolate solutes with elongation x > 1, although at more 
extreme solute volumes the values for prolate solutes are 
clearly lower than those for the corresponding oblate so- 
lutes. These variations are well described by the SPT 
prediction (fill) even though it was parametrized against 
spherical solutes only. Hence we conclude that the varia- 
tions in (3A/J.2 are mainly due to a combination of differ- 
ences in the mean radius of curvature i?2 , and the surface 
area 5*2. This is illustrated in Fig. |H1 which shows R2 and 
S2 as compared to the corresponding measures of spher- 
ical particles of the same volume. S2/R2 increases with 
increasing particle volume, and so eventually the vari- 
ation of /3A/i2 with x will mirror that of S2 shown in 
Fig.il 

At higher densities very similar trends are observed 
in the results from simulations and SPT. Results for 
p* = 0.5 and p* = 0.8 are shown in Figs. El and El 
respectively. At p* = 0.5 the simulation and SPT re- 
sults clearly show that a prolate solute has a lower excess 
chemical potential than an oblate solute with reciprocal 
elongation and the same molecular volume. This trend is 
not obvious in the simulation results at p* — 0.8 because 
we were only able to carry out calculations for V2 = V\ 
(although the SPT predictions for larger solute volumes 
will of course show this effect). Other trends apparent 
in Figs. [7| IH1 and ^| include the increase in (3Ap2 with 
molecular volume and, more significantly, the decrease 
in (5Ap2 with increasing p* or decreasing T*. Hence, 
the model system correctly shows that the solubilities of 
non-polar solutes in polar solvents (or at least those that 
can be described faithfully with a DHS model) will de- 
crease with increasing temperature. This variation im- 
plies that din K/dT* < 0, and hence that the energy 
change Au 2 < ©. Since d{T*\nK)/dT* < 0, it is 
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FIG. 7: Excess chemical potentials of hard-spheroid solutes 
with elongation x in DHS solvents with p* = 0.2 and various 
dipole moments: (left) p* = 0; (middle) p* = 1; (right) p* = 
2. In each case the symbols correspond to solute molecular 
volumes of Va = Vi (circles), V2 = 2Vi (squares), V2 = 5Vi 
(diamonds), V2 = WVi (up triangles), and V2 = 20 Vi (down 
triangles). The symbols are from MC simulations and the 
curves are the predictions of I1H . 
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FIG. 8: Fundamental measures of spheroids compared to 
the corresponding measures of spheres with equal volume, as 
functions of the elongation x: the mean radius of curvature, 
R2 (solid line); the surface area, S2 (dashed line). 
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FIG. 9: Excess chemical potentials of hard-spheroid solutes 
with elongation x in DHS solvents with p* = 0.5 and various 
dipole moments: (left) p* = 0; (middle) p* = 1; (right) ^i* = 
2. In each case the symbols correspond to solute molecular 
volumes of V% = Vi (circles), V2 = 2Vi (squares), V2 = 3Vi 
(diamonds), and V2 = 4Vi (up triangles). The symbols are 
from MC simulations and the curves are the predictions of 
itTTt . 

also clear that the entropy change AS2 < 10 • The 
magnitudes of the variations in (3Afi2 with temperature 
- and hence |Ai*2| and |As2| ~ increase with increasing 
solute molecular volume. These results are analogous 
to those for hydrophobic solvation, in which the pres- 
ence of a solute molecule in the solvent causes some sort 
of ordering of nearby solvent molecules with an associ- 
ated decrease of the solvent-solvent interaction energy. 
It should be emphasized that in the present case, there 
are no solute-solvent interactions beyond the hard-core 
repulsion, and so an energy change upon inserting a so- 
lute molecule in to the DHS solvent would arise exclu- 
sively from solvent-solvent interactions. In contrast to 
real hydrophobic solvation, however, we anticipate that 
there will be no temperature minimum in the solubility; 
the hard-core repulsions between model solutes and sol- 
vents are never small compared to the thermal energy, 
and so /3A/i2 will level off asymptotically to the values 
presented here for the systems with pf = 0. 



IV. CONCLUSIONS 

We have investigated the effects of molecular shape on 
the solvation of hard-core solutes in model polar solvents. 
The solutes are hard spheroids with a range of elonga- 
tions and molecular volumes, while the solvent consists of 
a fluid of dipolar hard spheres. Although these are very 
crude representations, the properties of the model sys- 
tems should bear some qualitative resemblance to those 
of real solutions. 

Using computer simulations we have calculated the dis- 
tributions of transient spheroidal cavities with specific 
elongations in the pure model solvents. For a given sol- 



FIG. 10: Excess chemical potentials of hard-spheroid solutes 
with elongation x and molecular volume V2 = Vi in DHS 
solvents with p* = 0.8 and various dipole moments: (left) 
p* — 0; (middle) p* — 1; (right) p* = 2. The symbols are 
from MC simulations and the curves are the predictions of 

G3- 



vent molecular dipole moment, the most probable cav- 
ity dimensions decrease with increasing solvent density 
reflecting the reduction in free volume. More interest- 
ingly, for oblate spheroidal cavities the shape correspond- 
ing to the maximum average radial semi-axis become less 
anisotropic with increasing solvent density; the same is 
true for prolate spheroidal cavities with respect to the 
maximum average polar semi-axis. This suggests that 
at high solvent densities, the probability of developing a 
non-spherical cavity is reduced due to tight packing of 
the solvent molecules. For a given solvent density and 
cavity elongation, the average cavity size decreases with 
increasing solvent dipole moment (or decreasing tempera- 
ture) due to the association of solvent particles and hence 
a reduction in the nearest-neighbor separation. 

The excess chemical potentials of hard-spheroid so- 
lutes of various molecular volumes have been calculated; 
these are related to the cavity distributions since there is 
no solute-solvent interaction beyond the hard-core repul- 
sion. For given solute molecular volume, solvent density, 
and temperature, the excess chemical potential is min- 
imized for spherical solutes, while this function is less 
for a prolate solute with elongation x than it is for an 
oblate solute with elongation 1/x. These observations 
are easily explained with reference to a simplified scaled- 
particle theory written in terms of the equation of state 
of the pure solvent, the solute-solvent interfacial tension, 
and fundamental measures of the solute molecule, i.e., 
the mean radius of curvature, the surface area, and the 
molecular volume. Firstly, the unfavorable interfacial 
tension contribution is minimized for spherical solutes. 
Secondly, the violation of x <-> l/x symmetry for a given 
molecular volume is due to the differences in the mean 
radius of curvature and surface areas between the two 
shapes. Since the solubility of a spheroidal particle is 
simply related to the excess chemical potential, we can 
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predict that for a given molecular volume, rod-like parti- 
cles should be more soluble than plate-like particles, and 
spherical particles should be the most soluble. To illus- 
trate the magnitude of this effect, consider a DHS solvent 
with p* — 0.8 and p* = 1, which roughly corresponds to 
a typical polar molecular liquid. For prolate and oblate 
solutes with moderate elongations of 3 and 1/3, respec- 
tively, and with molecular volumes equal to that of a 
single solvent molecule, Eqs. Q and ljTT|) with the values 
given in Table ITD predict that the solubility of the prolate 
solute will be about 17% higher than that of the oblate 
solute. Doubling the solute volume would result in the 
prolate molecule being approximately 35% more soluble 
than the oblate one. 

For given solute elongation and solvent density, the ex- 
cess chemical potentials of the solutes increase with in- 
creasing molecular volume, and decrease with increasing 
solvent polarity (or decreasing temperature). The for- 
mer trend is obvious, while the latter trend is due to the 



concentration of free volume mentioned above. The de- 
crease in solubility with increasing temperature is similar 
to that observed in the solvation of hydrophobic solutes 
in water, and reflects the fact that the entropy of solva- 
tion is negative, i.e., the presence of the solute molecule 
would lead to an ordering of the solvent. As a result of 
this enhanced ordering, the energy change is also neg- 
ative. In the present calculations, the negative energy 
change is entirely due to solvent-solvent interactions. 

Unfortunately, a rigorous test of the simulation and 
theoretical results against experimental data is not yet 
feasible due to the restricted set of available molecular ge- 
ometries and sizes; the key results presented here concern 
hard solutes of equal volume and with reciprocal aspect 
ratios. It would therefore be interesting to examine ex- 
perimentally the solubilities of, for example, nanoscopic 
colloidal particles with well controlled geometries and 
particle volumes that are large compared to the solvent 
so that the prolate-oblate asymmetry is enhanced. 
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